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We show that ion masses in superfluid 3 He ought to be 
enormously enhanced (by a factor of 10 2 ) as compared with 
the same ion masses in 4 He measured at low temperature. 
We calculate precisely the dependence of the effective mass 
on pressure in 3 He-B, and show that the coherent (ballistic) 
motion of ions in 3 He-B can be studied experimentally at 
T < (0.3 - 0.2) T c . 

PACS numbers: 67., 66.20+d 



The problem of ion motion in normal 3 He-liquid has 
been of long-standing interest, partly because of its con- 
nection with the " orthogonality catastrophe" , but mostly 
because theorists have had a hard time explaining it. 
The ion motion is greatly overdamped at low tempera- 
ture, by multiple scattering of 3 He quasiparticles, so theo- 
rists have concentrated on calculating the experimentally 
measurable ion mobility. Early perturbative calculations 
predicted a mobility fi(T) diverging as 1/T 2 below a 
temperature To = pp/M , where M is the bare ion effec- 
tive mass ( M ~ 100 — 260m3, depending on pressure, for 
the negative ion; here 1113 is the 3 He atomic mass). Ex- 
periments on both positive and negative ions 
flatly contradicted this prediction; /i_ is roughly constant 
through and below To, all the way down to the superfluid 
transition T c . 

However this problem is a strong-coupling one. The di- 
mensionless ion- 3 He coupling is g = ppat r /3ir 2 , with <j tr 
the transport cross section, and p > 1. The high-T scat- 
tering rate equals r = To g 3> To, which is why the ion 
motion is overdamped already for T>To. Moreover, it 
was realised by Josephson and Lekner P] that for T < T 
the ion recoil is not free, but Brownian diffusive, down to 
the unobservably low temperature T co h — To ge~ 9 . This 
diffusive motion means that it is meaningless to define 
an effective mass for the ion above T co h- The theory of 
ion mobility in normal 3 He has nevertheless been consid- 
erably refined since then [^0|-|l2|. 

One obvious way for experimentalists to see coherent 
motion of an ion in 3 He is to go to the superfluid phases, 
where the gap cuts off the " orthogonality catastrophe" . 
Remarkably, this possibility has not been explored, either 
in theory or experiment (although some mobility experi- 
ments have been done we return to these below). In 
this paper I give a detailed theory of ion dynamics, which 
is exact in the large-g limit. A very striking prediction 



emerges from this analysis - that the effective mass of 
ions in the superfluid phases will be very large (up to 
2 • 10 4 m3, or some 100 times the bare ion mass). I calcu- 
late the effective mass M e ^(P) as a function of pressure 
in the low-T limit in 3 He-B, and suggest how this pre- 
diction might be verified experimentally. This prediction 
(which is clearly out of the framework of the standard 
models p3| ) should constitute a very stringent test of 
our ideas of particle dynamics in a Fermi liquid. 

The Hamiltonian is that of a spherical object in a Fermi 
liquid environment: 



H = -MR 2 + H F + V , 
V = J dr V(r - R)p(r) , 



(1) 



where R is the ion coordinate, Hp the Hamiltonian of 
3 He, and p(r) the 3 He density operator. We make use of 
the path integral technique and integrate out the Fermion 
degrees of freedom fll],|l^] , and start by considering the 
case of normal 3 He. Using Feynman's path integral over 
R in imaginary time |I3] the effective action in the par- 
tition function can be written 
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drdr — J r n (R r — R r /)e 
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where Sq — J drA/R 2 /2, uo n — 2im/ (3 are Matsubara fre- 
quencies, and the influence functional, J- n {R), is related 
to the overlap integral between the initial and final Fermi 
liquid states with different local potentials [Of; one has 
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16tt 
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Here Sf — 5(R) and Si = 5(0) are the scattering 5- 
matrices at the Fermi energy in the final (the particle at 
the point R) and initial ( R = 0) states. The connection 
with the overlap integral (/ | i) is established by 



In I (/ I i) I 



dwcoth(w/2T) n 



TmT K (w,R);. 



(4) 



The effective action (^]3|) is correct provided we deal with 
heavy particles, M 3> TO3. 
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The formal expression (Q) is highly nonlinear in R and 
can not be solved in general. However in the strong cou- 
pling limit, g 1, we can restrict ourselves to a quadratic 
expansion 



which results in a simple quadratic action 

s w = MPj2("l + r\u n \)\R n 



(5) 
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where R T = ^3 Rne lw " T . Moreover, if we calculate 
the mean square value of the particle displacement us- 
ingEq.(§) 



P^((R T -R T 2 ) 
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we find |L2| that the expansion (j^) is justified in the nor- 
mal state down to T co h] the higher order terms in the 
expansion (pfR) 2 — Ca{pfR) a + ■ ■■ give rise to small 
corrections proportional to (20 /ng)C 4 In £ which can be 
neglected at T ^ T co h- The case of negative ions is of 
most importance here because for the hard sphere poten- 
tial with PfR- ^ 1 (R- is the bubble radius) the coeffi- 
cient C4 turns out to be very small, C4 ~ 10~ 2 — 10~ 3 . 

Now in the normal state , the mobility /i is given in 
linear response, and in the i? 2 -approximation, by 



n/e 



Muj 2 - f(u) M -iuj + T 



(7) 



where f(uj) — AitF r (uj) / R 2 , and e is the particle charge. 
While this describes diffusive motion for a; < T, one 
can also think of Eqs.(|6[0) as describing a frequency- 
dependent mass renormalization M e ^ /M = 1 + T/ \ 
uj n \. In the diffusive regime one assumes that M e f f < 
v 2 >oc T; then Einstein's relation a oc< v 2 > t/T, with 
r the scattering time, plus the experimentally observed 
/i(T) = const, leads to the conclusion that M e ^ oc 1/T, 
since certainly r > 1/T. This however is only very indi- 
rect evidence for a temperature-dependent effective mass. 

We therefore consider the ion motion in the superfluid 
state . To do this we must consider the overlap integral 
We start from the expression derived by Yamada 
and Yosida |l5|] at T — for the normal state which can 
be readily generalized for the case of finite temperatures 
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Here n w is the Fermi distribution function, G(oj) is the 
Green function in the initial state, and /SV is the scat- 
tering potential which distinguishes between the initial 
and final Hamiltonians, Ay = Vf — V{. We employ 
the standard definition of the retarded and advanced 
Green functions and their on-shell and off-shell parts: 
G R A = G* f ±G* n ; G = Gf fJ +tanh(o;/2T)G« . Note, 
that in the general case Green functions are matrices 
not only in the momentum space, but in the spin and 
electron-hole channels as well. 

Now we make use of the fact that the particle dis- 
placement is small compared to 1/pf- This means we 
may treat AV as a weak perturbation for arbitrary V: 

AVp p , = i(p - p')R (tOpp/ • (9) 
The leading i? 2 -term in Eq.(||) is given by the formula 
— I f f dEdE'{l-n E ,)n E 
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which can be easily rewritten in the form correspond- 
ing to Eq.(Q). Thus we can write the expression for the 
function f(u>) (see Eq. (]?])) as 

lmf u = np F dEdE'a(E,E')(n E '-n E ) S(E-E'-oj); 



a(E,E') = 
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where n = p F /3ir 2 is the particle density of 3 He. The 
static limit of Eq.dTJI) gives the mobility in the elastic 
model e/u = nppjdEa(E,E) {—dn E /dE) (see Ref. 

11). _ 

Untill now we have not specified the superfluid phase 
of 3 He, and Eqs.© are valid both for 3 He-A and 3 He- 
B. To observe the ballistic motion of ions experimentally 
we need the scattering time to be very long (r > 1 — 
10/is). This condition under any reasonable experimental 
arrangements may be satisfied only in 3 He-B at T <C T c 
(in fact T < (0.3 - 0.2)T C ). So, in the rest of this letter 
we concentrate on the effective mass calculation in 3 He-B 
at T = Q. 

It follows from the form of the effective action (^) that 
the mass renormalization is defined by the w 2 -term in 
the small frequency expansion of the functional integral 



SM = — 



dujXmf{uj) 
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(in 3 He-B the effective mass is isotropic). For T <^T c we 
can neglect the contribution due to normal excitations 
and further simplify Eqs.([ic|,p^]) to 
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where a(x, x') = a(x, x')/ttR 2 _ depends only on (ppR-), 
and the dimensionless frequencies x — E/A. In the sim- 
plest case of weak scattering potential (the strong cou- 
pling limit j > 1 still may be realized through a large 
number of weak scattering channels contributing to a) 
one can substitute the Green function in Eq.([To|) by its 
unperturbed value, which is equivalent to performing a 
u — v Bogoliubov transformation on the normal state am- 
plitude. After straightforward algebra we find 



6M U - V = — (1 + (72 2a tr ) — r — 
lb A 



(13) 



where a 2 = f do (1 — cos 2 9) da /do. This mass renormal- 
ization could be as large as 4 • 10 7713 at zero pressure. 
However, we demonstrate below that the exact calcula- 
tion for the hard sphere potential gives a value of 5M 
substantially different from Eq.fllg]). 

As pointed out in Ref. |r| the scattering matrix has a 
resonant behavior at energies near the gap edge which has 
to be treated exactly. First, we express conventionally 
the Green function in terms of the scattering T-matrix 
as G(u>) — V^T^yV^ 1 — V~ x and present the trace in 
Eq.(|i~0|) in the form 



Tr{T^ n {x)AV- l T^{x')AV- x } 



(14) 



where (AV -1 ) PP ' = i(p - p')R (V -1 )pp'. The analytic 
solution for the T-matrix was found in Ref. |L6]]. Since 
the energy spectrum of 3 He-B is spherically symmetric 
and does not depend on spin, and ap (where a is the 
fermion spin operator) is invariant under simultaneous 
rotations of the spin and momentum, it is clear that T- 
matrix is diagonal in the total angular momentum j and 
its projection m. Introducing the angular momentum 
eigenstates | j, to, l=j ± 1/2) =| jto±) one can represent 
the T-matrix as HJ13 



itN{0)T = T = 



h(K) -a 2 t 3 (-K)<r 2 
h(K) a 2 t x {-K)a 2 

*i = EE l s >< s l 4 s ; 

jm s— ± 

jm s — ± 



(15) 



-s >< s t{ s , 



where N(0) is the density of states in the normal phase, 
and t\ s and t J 3 s are known functions of frequency and 
phase shifts at the Fermi surface, Kj± = t&n5i=j±x/ 2 



Kj s (l ~ ipKj- a )/dj 



*3 S = (j>M KjsK^/djs 



(16) 



(1 + i P K js )(l - ipKjs) - (p/x) 2 K JS K^ s 



Here p(x) = x/(x 2 — l) 1 ^ 2 . Obviously, we have the same 
matrix structure for the on-shell retarded T-matrix as 
that in Eq. (|l5|) with the scattering amplitudes being re- 
placed by the on-shell ones t — > r 



-i0(\ x 
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(17) 



It is easy to show that the nonzero contribution to the 
trace comes from terms having the scattering potential 
only in combinations l/Vi-l/V v = wN(0)(l/ 'Ki~l/ 'K v ). 
So, we may conveniently replace the scattering potential 
in Eq.(|l4|) with the K-matrix. The final expression for 
cr(x, x') reads 



cr(x, x ) 



-1 



nppirR 2 



Tr{T R n (x)AK-%*(x')AK-i}, 



(18) 



There is one extra contribution to the effective mass 
due to bound states in the gap at Ej S = ±A cos(<5j+ — 
6j-) |16|. It is proportional to the number of occupied 
bound states multiplied by the 3 He quasiparticle mass 



SM h 



(19) 



The energies of the bound states approach the gap edge 
as A — Ej ~ A(d5j /dj) 2 /2. It is physically clear that 
states very close to the gap edge will disappear with the 
ion recoil in the scattering processes being taken into ac- 
count. Thus the maximum orbital number contributing 



to Eq.([19|) is defined by A — Ej ~ A/g (the recoil energy 
was calculated with the renormalized mass; see below). 
For the hard sphere potential with pfR- ^> 1 the phase 
shifts drop abruptly for j > ppR^ 1 which allows us to 
estimate the contribution of the bound states as 



5Mb =~ (pfR-) 2 TO3 ~ 37t<?to3 



(20) 



This mass is likely to be larger than the bare ion mass 
M, but still is much smaller than the renormalization 
defined by virtual transition in Eq.(|12|). 

The procedure of evaluating M e ff is straightforward 
now, because the trace determining the function <r(E, E') 
can be expressed entirely in terms of phase shifts at 
the Fermi surface which for the hard sphere potential 
are defined as tan Si — ii{pFR-)/ni(j>FR-), where ji 
and ni are the spherical Bessel and Neumann functions 
of order I. For any given pressure the set of parame- 
ters pf(P), A(P), and R-(P) allows to get the effective 
mass M e ff = M + 5Mb + 5M by numerical evaluation of 
Eqs.( p^| , ^8| ). In our calculations we used the normal state 
parameters taken from Wheatley's tabulation M|. The 
ion radius was taken from Ref. |T(| . Unfortunately, we 
found no tabulation for A(T — > 0, P) in 3 He-B, and had 
to rely on a weak coupling relation A(P) — a - 1.76 -T C (P) 
with the pressure- independent coefficient a = 1.12 ||. 
We think that A(P) is the most uncertain parameter in 
the present calculation. The bare ion mass is also un- 
known, but it is unlikely to contribute more than 10% to 
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the effective mass, and we simpiy neglected this contri- 
bution. 

Figure 1 shows the effective mass of the negative ion 
as a function of pressure. It was found to be as large as 
~ 2 • 10 4 TO3 at zero pressure dropping down to 3 — 4 • 



10 3 to 3 for P > 10 bar. The pcrturbative result, Eq.(13) 



is shown by the dashed line demonstrating the difference 
between the exact scattering amplitudes in 3 He-B and 
those obtained by applying the u — v transformation on 
the normal state amplitudes. 

The prediction of a very large effective mass should 
clearly be tested experimentally. There is already cir- 
cumstantial evidence for a large mass in the mobil- 
ity experiments of Nummila et al. ||, who found that 
elastic scattering theory seemed to explain their data 
down to the lowest temperature obtained in 3 Hc-B; this 
would be hard to understand using a bare mass assump- 
tion, since it would give a recoil energy 3> T already 
at T = OA T c . However the recoil energy calculated 
using the renormalized mass is much smaller, and an 
elastic scattering model is valid down to a temperature 
T el ^p%/(MT/A) = A/g. 

However what we really require is a direct experimen- 
tal test. One could search for resonant transitions be- 
tween the ion energy levels near the liquid-vapor inter- 
face [|o). The distance between the ion and the sur- 
face is large as compared to the coherence length in 
3 Hc-B up to the electric field strength E ~ lOOV/cm, 
with a typical range of resonance frequencies around 
ljq ~ 10 — 40 MHz. From the mobility experiment || we 
estimate that loqt 3> 1 below 0.3 T c , and the resonance 
is sharp enough to be observed. In the time-of-flight ex- 
periment the ion mass in the bulk can be measured for 
arbitrary pressure. In this technique the electric field 
is reversed during a time interval At ~ 1 — 10/is, and 
the ion is supposed not to accelerate above the Landau 
critical velocity eEAt/M e ^ < A/pf- Fortunately, the 
scattering time is as long as 1/is already at 0.3T C , and the 
enormous ion mass effectively compensates the smallness 
of the critical velocity in 3 Hc-B. 
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Effective mass of the negative ion in 3 He-B. The solid 
line is the exact calculation for the hard-sphere poten- 
tial which is compared to the perturbative result (dashed 
line) with the same transport cross section in the normal 
state. 
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